SCIENTIFIC 

REPORTS 




^ Acoustic rainbow trapping 



SUBJECT AREAS: 
APPLIED PHYSICS 
MECHANICAL ENGINEERING 
FLUID DYNAMICS 
MECHANICAL PROPERTIES 



Jie Zhu 1 , Yongyao Chen 2 , Xuefeng Zhu 1,3 , Francisco J. Garcia-Vidal 4 , Xiaobo Yin 1 , Weili Zhang 2 
&Xiang Zhang 1 

] NSF Nano-scale Science and Engineering Center (NSEC), 3112 Etcheverry Hall, University of California, Berkeley, California 
94720, USA, 2 School of Electrical and Computer Engineering, Oklahoma State University, Stillwater, Oklahoma 74078, USA, 
3 Huazhong University of Science and Technology, Wuhan, Hubei 430074, China, 4 Departamento de Fisica Teorica de la Materia 
Condensada and Condensed Matter Physics Center (IFIMAC), Universidad Autonoma de Madrid, 28049 Madrid, Spain. 



Received 
21 January 2013 

Accepted 
4 April 201 3 

Published 
25 April 201 3 



Correspondence and 
requests for materials 
should be addressed to 
W.L.Z. (weili.zhang@ 
okstate.edu) or X.Z. 
(xiang@berkeley.edu) 



Spatial modulation of sound velocity below the wavelength scale can introduce strong frequency-dependent 
acoustic responses in tailored composite materials, regardless the fact that most natural bulk materials have 
negligible acoustic dispersions. Here, for the first time, we experimentally demonstrate a metamaterial that 
traps broadband acoustic waves and spatially separates different frequency components, as the result of 
dispersion and wave velocity control by designed gradient subwavelength structures. The trapping positions 
can be predicted by the microscopic picture of balanced interplay between the acoustic resonance inside 
individual apertures and the mutual coupling among them. With the enhanced wave- structure interactions 
and the tailored frequency responses, such metamaterial allows precise spatial-spectral control of acoustic 
waves and opens new venue for high performance acoustic wave sensing, filtering, and nondestructive 
metrology. 



Slowing down optical waves with resonating photonic structures introduces controllable optical delays and 
allows temporary storages of light 1 " 5 , resulting strong light-matter interactions lead to novel optical devices 
such as all-optical memories and switches. The causality, however, fundamentally restricts the attainable 
delay-bandwidth product, limiting the resonance-based slow- light applications to a narrow band of frequency. 
Adiabatic control of the dispersions lifts such a constraint and allows the recent demonstration of slow light 
device 6 ' 7 and optical trapped rainbow 8 " 15 . In stark contrast to the tremendous progress in slow-light photonics, 
decelerating, trapping and spectrum-splitting of acoustic waves have not yet been realized due to the fact that 
nature occurring materials lack strong acoustic dispersion, which however is critical in many applications ranging 
from cochlear bionics to frequency mixing and selections in acoustic wave detection and imaging. Much like their 
optical counterparts such as photonic band-gap crystals and metamaterials 16 " 18 , artificial acoustic materials have 
shown promising potential of manipulating the acoustic wave and the wave-material interactions below the 
wavelength scale 19 " 21 . Innovative ways of acoustic wave collimation, focusing, cloaking, imaging, and extraord- 
inary transmission have been explored 22 " 27 . However, strong enough acoustic dispersion that slows down and 
eventually leads to the trapping of sound waves with artificial materials has only been theoretically proposed 28-30 . 
By utilizing metamaterial design and tailoring its dispersion to an extreme level, here we experimentally dem- 
onstrate a new class of anisotropic metamaterials that can efficiently trap broadband acoustic waves and spatially 
split different frequency components. This effect is due to strong modulation of wave velocity through gradient 
subwavelength unit cells that are strongly coupled along the propagation direction. The decelerating process can 
be theoretically explained within the effective medium approach, which however tends to break at the trapping 
point. The accurate trapping locations are otherwise predicted by interpreting a balanced microscopic interplay 
solution, as we present in the following. 



Results 

The acoustic rainbow trapping metamaterials, as shown in Fig. la, consists of an array of grooves perforated on a 
rigid bar. The grooves, with linearly increasing depth h, have width w and period p that are much smaller than the 
wavelength of interest. This 2D material is described as an effective medium model of a tapered anisotropic 
metamaterial layer (medium I) sandwiched between the homogeneous and isotropic layers of environmental 
medium II and rigid body III (Fig. lb). When considering an acoustic waves propagating along the metamaterial, 
we can analytically derive the corresponding group velocity v g of the fundamental guided mode from the general 
dispersion equation (see Supplementary Equation S8 online): 
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Figure 1 | Metamaterial for acoustic rainbow trapping, (a) The 

metamaterial is made of 80 grooves perforated on a square brass alloy bar, 
with period p = 6.35 mm and width w = 4.76 mm. 80 grooves are equally 
divided into 16 groups, which each contains 5 uniform grooves, while h of 
each group increases from 1.59 mm to 25.4 mm with step of 1.59 mm. 
(b) Effective medium model of the composite structure consists of 
anisotropic metamaterial medium I sandwiched between two 
homogeneous and isotropic media II and III (respectively air and rigid 
material in this study). Orange arrow indicates wave propagation 
direction. In ideal two dimensional studies, the length of waveguide 
extends infinitely in the y direction. Due to the limitation of fabrication 
technique, actual metamaterial extends 44.45 mm, which is multiple times 
over largest wavelength of interest, in y direction. 
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Here k 0 = w/c air is the free space wave vector. 0, as shown in Fig. 2b, 
is the tapering angle and z is the distance along the z axis. The 
effective anisotropic mass density and bulk modulus of metamaterial 
layer I are respectively p = (p/w)p ain p z = 00 and k = (p/w)K ail , As 
shown in Fig. 2a, calculated group velocity of the acoustic wave with 
given frequency /progressively slows down and eventually stops at 
z — c air /(40f), as implied by equation (1). Such spatial evolution of 
acoustic wave group velocity is hardly found in nature materials, 
revealing the metamaterial's strong acoustic dispersive behavior, 
which is evidently shown in Fig. 2b with significantly energy con- 
centration of three different spectral components (5 kHz, 7 kHz and 
9 kHz) at various locations. This result is validated by the full wave 
2D FEM simulation presented in Fig. 2c, suggesting that broadband 
incident acoustic waves can be trapped into spatially separated loca- 
tions and therefore spectrally split, resembling the phenomenon of 
optical "trapped rainbow" 8 " 15 . 

The experiments were conducted in air with a setup schematically 
shown in Fig. 3a. Figure 3b shows the position- dependent acoustic 
wave transmission along the metamaterial, which is measured at 4.5- 
9.5 kHz. Different frequency components are gradually stopped and 
successively localized along the waveguide, progressively inducing 
enhanced acoustic fields that approach the maximums. The intens- 
ities measured at the trapping point have over two orders of mag- 
nitude boost over non-trapping scenario, where the acoustic field is 
measured at the same distance from the speaker with reference metal 
bar containing uniform grooves. Such results are well agreed with the 
3D full wave FEM simulation shown in Supplementary Fig. S2 
online. The minor differences arise from the deviation that happens 
when simulation setup tries to accurately mimic the actual speaker 
configuration and ambient environment. The spatial evolutions of 
acoustic pressure fields for various frequencies (4.5 kHz, 6.0 kHz, 
and 9.5 kHz) are illustrated in Fig. 3c. Alike one can be concluded 
from the dispersion relation in Fig. 2a that is derived based on the 



effective medium model, higher frequency acoustic waves are clearly 
trapped in the region with shallower grooves, whereas lower fre- 
quency component like the 4.5 kHz one propagates along the meta- 
material until stops at section with deeper grooves. The experimental 
results agree well with full-wave 2D finite element method (FEM) 
simulation conclusion (Fig. 4b). It can be observed that the experi- 
mentally measured spectral splitting ( frequency / versus trapped 
position z ) exhibits a nonlinear distribution along the linearly 
tapered metamaterial, implying the non-uniform spatial density of 
the frequency components \df stop /dz\ =c air /(40z 2 ), which can be 
derived from equation (1) and noticed in Fig. 2a. However, by tailor- 
ing the position-dependent groove depth gradient, arbitrary spatial 
separation of various spectral components can be realized. 

To achieve the strong dispersion required for acoustic rainbow 
trapping, the groove fins of the metamaterial should be prevented 
to support any mode that could disturb the acoustic field in air. We 
have accordingly picked brass alloy as the base material, taking 
advantage of the modulus of rigidity about six orders of magnitude 
larger than that of air. Thus the influence of fin vibration is negligible 
when we experimentally study the fabricated metamaterial structure 
in solid-fluid system. We also conducted multiple measurements 
using the same setup. The results show small variations in acoustic 
pressure fields and trapping locations that result in the error bars in 
Fig. 3. Here it is worth noting that the demonstrated broadband 
acoustic wave decelerating, trapping, and spectral splitting do not 
rely on any phase manipulation along the propagation direction, very 
different from previously proposed acoustical cochlear-like response 
and place-to-frequency maps 31 " 33 . 

The effective medium model can well describe the propagation 
and gradually slowing down of acoustic waves because of the deep 
subwavelength nature of unit cells in the metamaterial. However, 
when approaching the trapping location, the effective wavelength 
of acoustic wave is progressively decreased. Therefore these unit cells 
are no longer deep subwavelength in size, which leads to the break- 
down of effective medium model, the deviation at predicting where 
the acoustic wave will be trapped and how the intensity of acoustic 
wave is distributed around the trapping location, as witnessed by 
different intensity peak positions and dropping rates observed in 
Fig. 2 and Fig 3. Under such circumstance, the mechanism of the 
acoustic wave localization at trapping point can be better interpreted 
by a microscopic model that the acoustic wave trapping originates 
from the balanced interplay between the acoustical oscillation inside 
individual grooves and the mutual near-field coupling among the 
neighboring grooves, as depicted in Fig. 4a. This coupling is governed 
by the diffracted evanescent waves with maximum momenta ± nip 
that propagate along the interface between the metamaterial and the 
air, driving the pressures and particle velocity fields to oscillate in 
opposite phases in neighboring grooves, leading to a standing wave 
with diminishing group velocity at the trap. As a result, the accu- 
mulation of the acoustic energy and field enhancement is observed in 
Fig. 3. A closed-form-expression for this model can be derived (see 
supplementary information for details): 



cot (k 0 h)-- 



Min(ff) 



(2) 



Here cot(k 0 h) defines the acoustical oscillation localized within 
individual grooves, while the right side term is dominated by the 
diffracted evanescent waves, relating to the coupling effect discussed 
above. Balance between those two physical processes results in the 
trapped wave. Microscopically in Fig. 4a, it can also be seen as a 
series of closed-loop circulations of acoustic energy-flux between the 
adjacent grooves, resulting in the standstill of acoustic wave with 
zero group velocity. It's worth noting that for this microscopic 
model, only under the strict condition when the period of unit cells 
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Figure 2 | Spatial modulation of acoustic waves investigated within the effective medium approach and full wave simulation, (a) Calculated group 
velocities (normalized to the speed in air) of acoustic waves (5 kHz, 6 kHz, 7 kHz, 8 kHz and 9 kHz, respectively) propagating along the metamaterials 
shown in Figure 1. Various spectral components are progressively slowed down and eventually trapped at different locations which are marked with 
colored arrows. The spatial density of trapped spectral components (space between the colored arrows) is non-uniform along the linearly tapered 
metamaterial. (b) Simulated spatial energy distribution of different acoustic waves (5 kHz, 7 kHz, and 9 kHz), presented in RGB colors (Red, Green and 
Blue), shows significantly enhanced acoustic fields at the positions where different spectral components are trapped. 6 is the tapering angle of 
metamaterials core. Simulation was carried out with partial differential equation module. The wave intensity is presented by the strength of each 
individual color, (c) Acoustic field intensity variation from full wave simulation on 2D actual structure reveals the same tendency of gradual field 
enhancement as the one derived from effective medium approach. 



p — > 0, there is an asymptote 3tf sto p = c air /4h derived from equation 
(2), and each groove can be treated as a quarter-wave resonator for 
trapped acoustic wave, which is consistent with equation (1) as per- 
fect effective medium. However, due to the finite size of actual unit 



cells, the microscopic model is more properly suited to tell the exact 
acoustic wave trapping position and agrees better with the 2D full- 
wave FEM simulations on real metamaterial geometry than effective 
medium model, as presented in Fig. 4b. 
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Figure 3 | Experimental demonstration, (a) Schematic drawing of the experimental set-up. A guiding metal bar (dark Gray) is introduced to excite 
bound modes propagating along the surface, coupling acoustic waves towards the metamaterial sample, (b) Normalized acoustic fields intensity versus 
the measurement positions (distance from the start of metamaterial in z direction) shows clear rainbow trapping effect, with incident waves changing 
from 4.5 kHz to 9.5 kHz by step of 25 Hz. (c) Acoustic field intensity variation along the metamaterial for 4. 5 kHz, 6.0 kHz, and 9.5 kHz incident waves. 
All three show progressively enhanced fields before peaking at spatially separated locations. The error bars represent s.d. among 5 repeat measurements. 



Discussion 

In summary, we have experimentally demonstrated that broad- 
band acoustic waves can be gradually slowed down to standstill, 
thus trapped and spatial-spectrally separated when propagating 
along a metamaterial waveguide. Such "trapped rainbow" effect 
is achieved by the dispersive modulation of acoustic velocity 
through gradually varying subwavelength unit cells that are 
strongly coupled along the propagation direction. The acoustic 
"trapped rainbow" can be released by disrupting the balance 



between the acoustical oscillation inside individual grooves and 
the mutual near-field coupling among the neighboring grooves. A 
simple example is to place a small object at the top or inside 
trapping groove to strongly scatter otherwise trapped acoustic 
waves. Contrary to the optical counterparts 34 , this acoustic meta- 
material with the unique spectral splitting capability is also suit- 
able for 3D system due to the absence of cutoff wavelength. 
Without bandwidth limitation, it opens the door to realizing pre- 
cise spatial-spectral control of acoustic waves. 
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Figure 4 | Microscopic mechanism associated with steady acoustic wave trapping state, (a) Acoustic wave trapping is the result of balanced interplay 
between the local groove oscillation and the mutual near field coupling among the grooves to form closed-loop circulations of "trapped sound". When 
acoustic waves are trapped, the effective acoustic energy flow (black arrows) between two adjacent grooves will go opposite directions during 0-t/2 (top 
half) and t/2-t (bottom half), here t is the period of oscillation inside individual groove, thus forming an equivalent closed-loop circulation where the 
acoustic wave always returns to its origin, resulting in the standstill of acoustic waves. The vertical white arrows represent acoustical oscillation inside the 
grooves. The horizontal white arrows correspond to the diffracted evanescent waves that propagate along the interface between the metamaterial and the 
air. The color map is the normalized acoustic energy flow, (b) Comparison of acoustic wave trapping locations derived by 2D full wave FEM simulation on 
actual metamaterial geometry, the microscopic model, the quarter-wavelength resonator model, and experimental results. Microscopic model result fits 
well with the FEM simulation curve and experimental results, suggesting better capability in analyzing the steady acoustic wave trapping state, especially 
when the unit cell size is not infinite small. Experimental results start to deviate from theoretical prediction at lower frequency range, which is due to the 
larger wavelength and finite length of actual sample in y direction. The error bars represent s.d. among 5 repeated measurements. 



Methods 

Sample preparation. We fabricated the metamaterial structure with a square brass 
alloy bar. 80 grooves are perforated on the bar by high pressure water jet cutting, with 
periodp= 6.35 mm and width w = 4.76 mm. 80 grooves are equally divided into 16 
groups, which each contains 5 uniform grooves, while h of each group increases from 
1.59 mm to 25.4 mm with step of 1.59 mm. 

Experiments setup and measurements. Experiments are conducted under room 
temperature inside a large chamber to isolate the environmental noise. Attached on 
all the interior walls, cone shaped sound absorbing foam resembles the perfect 
matching boundary around the setup. A square brass alloy bar with same width as the 
metamaterial and uniform shallow grooves, connected with the metamaterial at one 
side and acting as an ordinary input waveguide, couples the incident acoustic waves 
from the speaker into the metamaterial. A 17 mm diameter PUI Audio AST- 
1732MR-R round speaker is located at the other side of input waveguide, sending out 
continuous sinusoidal waves towards the sample at various frequencies. The diameter 
of speaker is smaller than the shortest wavelength applied in test so that various 
spectral outputs possess no directivity difference. Incident signal were generated by 
Tektronix AFG3021B Arbitrary/Function generator. Acoustic pressure fields are 
measured by inserting a 4 mm diameter PUI Audio TOM- 1045S-C33-R microphone 



from the side into each of the 80 grooves and close to the bottom. This microphone is 
attached on a thin metal wire to minimize field perturbations. All measured data was 
sent to computer through Tektronix TDS2002B digital storage oscilloscope for 
further processing. Both AFG3021B and TDS2002B were controlled by Lab VIEW® 
software. 

Numerical calculations. FEM simulations on acoustic pressure field measurement 
were carried by COMSOL Multiphasics™ 4.1 with the acoustic- solid interaction 
multi-physics module. The materials applied in simulations were air and ultra- 
machinable brass (Alloy 360). Plane wave radiation boundary conditions were given 
to the outer boundaries of simulation domain so there will be no interference from 
reflected acoustic wave. The largest mesh element size was set lower than 1/10 of the 
lowest wavelength, and finer meshes were applied at the domain of metal fins. 
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